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Field-induced magnetic order has been investigated in detail in the in- 
teracting spin 3/2 dimer system Cs3Cr2Br9. Elastic and inelastic neutron 
scattering measurements were performed up to = 6 T, well above the criti- 
cal field Hci ~ 1.5 T. The ordering displays incommensurabilities and a large 
hysteresis before a commensurate structure is reached. This structure is fully 
determined. Surprisingly, the lowest excitation branch never closes. Above 
Hci, the gap increases slowly with field. An analysis in terms of projected 
pseudo-spins is given. 

PACS numbers: 75.10.Jm, 75.40.Gb, 78.70.Nx 



In recent years, a great deal of work has been devoted to isotropic quantum spin systems 
which show a gap in their energy spectrum. This includes a large variety of different models 



such as Haldane and alternating chains, ladders, as well as spin-dimer systems. In zero 
magnetic field, these systems are all characterized by a 5" = singlet ground state (|0)) and 
the energy gap Eq of the lowest magnetic excitation corresponds to a 5* = 1 state In 
a field H, when the gap closes at if = Hd {qhbHci ~ Eq) a three-dimensional magnetic 
ordering develops, induced by the small inter-chain and/or inter-dimer couplings. In these 
field-induced magnetic orderings (FIMO) only the spin components transverse to the applied 
field spontaneously break symmetry. Moreover, as initially proposed for Haldane chains, a 
FIMO transition is an experimental realization of a Bose-Einstein condensation (BEG) of 
the hard-sphere bosons - the magnons - into the ground-state [1]. Recently, the magnetic 
behavior observed in the s — \ spin-dimer compound TICUCI3 has been analyzed within 
this framework [2,3]. While for BEG a gapless Goldstone mode is present in the energy 
spectrum of the condensed phase, in the case of a FIMO, the mode is gapless only if there 
is conservation of the component of the total spin in the field direction. If this symmetry 
is broken, the total number of the magnons is no longer conserved. A gap may then be 
seen in the lowest excitation branch. It is important to understand the way an anisotropy 
creates such a gap in a FIMO phase. In this letter, we present a complementary elastic 
and inelastic neutron investigation of the FIMO phase in the spin-dimer system Cs3Cr2Br9. 
This material differs from TlCuCls and similar compounds [4] in two significant respects. 
The larger spin value of Cr^+, s — 3/2 favors the presence of local anisotropics. The crystal 
structure is different: the spin-dimers form an hexagonal arrangement in the (a, b) plane as 
shown in Fig. 1. This leads to frustration, which should strongly affect any magnetic order. 
These two points, i.e., the role of a small anisotropy and the effects of frustration in a FIMO 
are, here, demonstrated experimentally. An analysis in terms of projected pseudo-spins is 
proposed to discuss the striking features in the magnetic phase. 
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FIG. 1. Transverse spin order of the Cr^"*" ions in the commensurate + H magnetic phase of 

Cs3Cr2Br9 for H || a — b. k is the propagation vector of the magnetic order. 

The spin-excitation spectrum of Cs3Cr2Br9 in zero field was determined across the whole 
Brillouin zone by Leuenberger et al. [5]. The intra- and inter-dimer exchange couplings are 
all antiferromagnetic (AF). In a field, Zeeman splitting was also observed, but only for small 
fields, H < Hci [6]. Neither the field- induced magnetic structure nor the field dependence 
of the spin excitations have been studied above Hd- Here, we investigate this FIMO phase 
in detail. Our experimental configuration is the same as in Ref. [6]: the external field H is 
applied within the (a, b) plane along the a - b axis (hereafter defined as the Z axis) so that 
the XY plane where the AF ordered moments are expected to lie {transverse ordering) is 
defined by the c (hereafter X) and the a -|- b (Y) axes (Fig. 1). Cs3Cr2Br9 crystallizes in 
the hexagonal P&^/mmc structure, with a = b = 7.508 A, c = 18.70 A (which is twice the 
distance between successive planes), and consists of four Bravais sublattices. The s — 3/2 
Cr^+ dimers lie parallel to the c axis [7] . 
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FIG. 2. a) Field-dependence of the propagation wave-vector component k {= kx = ky) and b) 
that of the parallel (Mz) and the ordered transverse (Mx) magnetizations. Qjj ^ is the wave 
vector corresponding to the minima of the dispersions (see Fig. 3 and text). In panel 2b, the solid 
lines are theoretical predictions obtained for the commensurate + H magnetic order (see text) . 

First, we present the elastic results. They were performed on the lifting-arm two-axis 
diffractometer D23 CEA-CRG at the Institut Laue-Langevin (ILL) in Grenoble (Prance). A 
single crystal (~ 350 mm^) of Cs3Cr2Br9 was placed in a 6T vertical cryomagnet equipped 
with a dilution insert. The wavelength was A = 1.28 A and the temperature was maintained 
in the range 50 — 100 mK. In zero field, 113 nuclear reflections were collected, which could 
be reduced to 68 independent measurements. The structural refinement was achieved with 
type II Becker-Coppens Gaussian extinction correction and isotropic thermal factors. A 
weighted R-factor of 1.7 % on the nuclear structure factors Fjv was obtained by refining 12 
parameters. The magnetic structure was then determined at if = 6 T. It is described by 
the commensurate propagation wave vector k = (1/4, 1/4, 0) [8]. Due to the symmetry, 
two other magnetic domains were observed, corresponding to the wave vectors (1/4, —1/2, 
0), ( — 1/2, 1/4, 0). 123 magnetic refiections were collected in the three domains. For the 
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magnetic refinement an isotropic Cr^+ magnetic form factor was used, yielding a weighted 
R-factor of 4.2 % on the magnetic structure factors Fm- This gives: i) the populations of the 
three domains are about the same (34.8 ± 0.6, 33.0 ± 0.4 and 32.2 ± 0.4%) and the magnetic 
structure is coUinear with the moments pointing along the c (= X) direction; this agrees with 
the expected transverse ordering for a FIMO, ii) the arrangement of the transverse magnetic 
moments within one cell is represented in Fig. 1 (atoms labelled 1, 2, 3, 4). If we assume the 
moment amplitude to be constant (as is usual for 3d ions), this motif propagates with the 

sequence + H along the (110) direction (for the first domain), as shown in Fig. 1, and 

after refinement, we obtain that the AF ordered transverse moment is M.x = 1.60(1) iib per 
Cr^+ ion at H — 6T. Considering the black atoms in this figure, ferromagnetic correlations 
between neighboring spins, i.e., frustration, are seen both within the (a,b) plane and between 
the planes. The field-dependence (for i7 < 6 T) of the propagation wave-vector component 
k — kx — ky and that of A4x are displayed in Fig. 2 (open and sohd circles). These data 
were from scans in the (110) direction across the strongest magnetic peak (1 — A;, 1 — k, 2) 
with k ~ 0.25, for increasing and decreasing magnetic field H. When increasing (decreasing) 
H, this peak could be detected from (down to) if ~ 1.5 T: Hd is, then, slightly smaller 
than this value. Each scan was well fitted by a Gaussian function, giving rise to an accurate 
evaluation of both the integrated intensity, i.e., /Ax^ the wave- vector component k (see 
the inset in Fig. 2a where the magnetic peak is recorded for the same field value, H — 2.9 
T, while H is increased or decreased). In Cs3Cr2Br9, within the ordered phase, there are 
both incommensurabilities and large hysteresis effects suggesting pinning [9]. For increasing 

H, the commensurate -|- H structure is reached at if ~ 5 T, while for decreasing H, 

the commensurate structure remains locked down to if ~ 3 T. The parallel magnetization 
A4z{H) was measured both by neutrons from the ferromagnetic contribution appearing on 
top of the weak nuclear Bragg peak (1, 1, 4) at if = 6 T (the open triangle in Fig. 2b) and 
from SQUID at T ~ 100 mK in the full field range (solid triangles). 

Second, we report inelastic results. The same crystal was mounted on the cold neutron 
triple axis spectrometer IN12 FZ/CEA-CRG at ILL, in a 6T vertical cryomagnet with a 
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dilution insert. The measurements were performed in the same experimental configuration 
at T ~ 50 mK. The final wave-vector \s.f was fixed to various values between 1.07 and 1.55 
A~^, depending on the explored energy range and on the energy resolution that was needed. 
Examples of the dispersion as a function of Qh,k — Qh — Qk (at Ql — 2) are shown in 
Fig. 3. 
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FIG. 3. Examples of dispersions for H = and H = 3 T. The solid line for if = is a fit as in 
Ref . [5] . The dashed lines for = 3 T are guides to the eye. 

For H — (open squares), a unique excitation branch is observed, while three distinct 
branches are seen in the ordered phase (sohd circles for if = 3 T) . The description proposed 
in Ref. [5] for if = reproduces very well the data (sec solid line). From this fit, we obtain 
the following exchange values: J = 1.03 meV for the intra-dimer coupling, Jp — 0.054 and 
Jc — 0.039 meV for the in-plane and out-of-plane inter-dimer interactions, respectively, in 
agreement with Ref. [5]. We now focus on the minima of the dispersions, i.e., the three 
energy gaps. They always occur at the same wave-vector value Qhk ~ 0.27 r.l.u., even 
above ifci where the propagation vector component k varies continuously from 0.29 down 
to 0.25 r.l.u. (see Fig. 2a). The field-dcpcndcncc of these gaps, obtained from energy scans 
performed at (0.27, 0.27, 2), is plotted in Fig. 4. The inset shows such a scan recorded at 
H — with the best energy resolution (kf — 1.07 A"-*^): a splitting - not detected previously 
[5] - is clearly observed [10]. It reveals the presence of an additional anisotropy which can 
be described by a single-ion term D[sfY (X = c) at each Cr^+ site. Such a term results in 
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a splitting of the initial dispersion. This effect can be accounted for by redefining for each 
branch two distinct effective intra-dimer couphngs [11]: J" J + 0.8D and = J — 1.6D 
(for s = 3/2 and D < 0). From the splitting at the minimum of the dispersions, we evaluate 
D ~ —0.01 meV. With this value, the Zeeman splitting below Hd (solid lines in Fig. 4) 
agrees also very well with the data. From the lowest energy branch, the critical field is 
expected at Hd ~ 1.5 T, in agreement with the diffraction result. As mentioned above, for 
H > Hci, three distinct branches are still observed. The two upper branches agree roughly 
with "Zeeman" behaviors: qhbH and ^qiibH with g = 2 (dotted hues). The evidence for 
a gap {£g) on the lowest energy branch above Hd is, however, a new result. It shows that 
there is no massless Goldstone mode in the FIMO phase of CszCr2Br^. In fact, even near 
Hd the lowest energy gap does not vanish. The half-width (0.05 meV) of the fitted gaussian 
at the minimum gap ~ 0.11 meV is instrumental. 




FIG. 4. Field dependence of the energy gaps. The solid lines describe the splitting below Hd- 
Above Hd, the dotted lines are gixsH and 2giJ,BH. The dashed line describes the lowest gap Eg 
(see text). Inset: H = energy scan performed at (0.27, 0.27, 2)with kf = 1.07 A~^. The splitting 
signals the presence of a small single-ion anisotropy {D in Eq. 1). The slightly broader width of 
the upper mode indicates the possible presence of small higher-order anisotropies. 

Including the single-ion anisotropy D, the Hamiltonian for the spin-dimers in Cs3Cr2Br9 

is 

n = Y.[J^ia-Sib + D{sf," + ) - 9^^BH{sl + 4)] 

i 
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~l~ ^ ] J ipia-^ma ~l~ ^ib-^mb) (l) 

with J' — Jp (Jc) for in (out-of)-plane couplings. In Eq. (1), the indices a and b refer to the 
two spins in each dimer while the index m accounts for the interacting spins surrounding 
each spin i (6 in-plane and 2x3 out-of-plane couplings). Because of the spin value s = 
3/2, three distinct plateaus are, a priori, expected in the magnetization curve of such a 
compound. They are actually observed in the parent compound Cs3Cr2Cl9, and magnetic 
ordered phases are expected to precede each of the three plateaus [12]. The situation is 
different in Cs3Cr2Br9 where the magnetization constantly increases from Hd ~ 1.3 T up to 
saturation {Hg ~ 26 T) without intermediate plateaux [12]. This means that the three FIMO 
phases in this compound overlap giving rise to a single FIMO phase, hereafter called the 
extended FIMO phase. In general, a simple FIMO phase develops each time a crossing occurs 
between two states of the isolated dimers. In low field, H ~ Hd, the crossing occurs between 
the singlet ground-state |0) and the lowest Zeeman split state |1, — 1) (including anisotropy) 
of the initial triplet state |1). For a simple FIMO phase, an effective Hamiltonian Hgff 
can be derived by projecting Eq. (1) onto these elementary states. One obtains & a — \ 
pseudo-spin hamiltonian, which, for the initial s = 3/2 spin system, is given to first order: 



n-l \ ' tXY (XX , „Y Y \ , jY „zz 

T^eff = 2^ Jeff {(^i + ^mJ + Jeff^i 

i,m 

+deff{crf(^^ - crfa^) - ^ g^ieHeffaf (2) 

i 

With J^^^.=2M^J', J^ff = J'{cl-c\f/2, deff^ -AM^J'cid and giieHeff = gi^BH{cl - 
di) + J+J'{cl-c\)/A-R^D. In these expressions, M= °(1, -l|sjf |0) = — ^5/2, R=-2{1 - 
6ciC_i)/5, Ci and c_i account for the mixing of the initial states (no anisotropy): |1, —1) = 

r n 1/2 

c_i|l, -ly + ci|l, 1)° with c_i = (3D I [{(3Df + {[{giiBHf + {(3DfY''' - gi^BHY] and 
ci = {[{g^BHfH(^DfY'^-9^^BH}/ [{^Df + {[{g^iBHf + {^DfY'^ - g^^BHY] where 
/3 = 2 ''(l,-l|s^|l,+l)'' = 6/5. We note that the ratio J^fflJ^f^f ~ 1/(4M2) = 1/10 is 
much smaller than in the s — 1/2 spin case [13], i.e., for s — 3/2, the pseudo-spin hamiltonian 
is almost XY. More novel is the XY anisotropy term d^ff- This anisotropy is induced by 



D, but it occurs indirectly through the mixing coefficient cic_i. Note that higher-order 
anisotropics within the YZ hexagonal plane (see Fig. 4 caption) would renormalize slightly 
the ratio Jeff/'^eff deff, but otherwise does not affect the analysis. According to 
Eq. (2), in a simple FIMO phase, one has an anisotropic quantum pseudo-spin 1/2 XYZ 
model. For an extended FIMO phase - as in Cs3Cr2Br9 - the pseudo-spin description should 
give the essential physics, but one expects deviations. We now compare our data and 

predictions from Eq. (2). For the a spin system, we analyze the commensurate + -\ 

structure of the FIMO phase. The effective couphngs in the mean-field approximation 
according to the spin arrangement displayed in Fig. 1 (see the black atoms) are {Jeff) — 
-4(Jp + Jc/2)M'^ and (4//) = +4.(Jp + Jj2)M'^ciC^i. The calculated Mz and Mx, with 
M.X — 2|M|y'A42(l — Mz)-, are represented by the solid lines in Fig. 2b. The observed 
behaviors are well explained qualitatively. Quantitatively, at i7 = 6 T, the discrepancy 
between measurements and this simple theory is ~ 15 % for Aix- Such a correction is what 
one expects when one compares the results for a simple and an extended FIMO, i.e., when 
the contributions from the higher states are taken into account [14]. Extrapolation of these 
curves to lower fields predicts a reasonable critical value, Hci ~ 1.2 T. 

A FIMO corresponds to an AF spin- flop state. A small anisotropy, such as dg// in Eq. 
(2), breaks the axial symmetry responsible for the Goldstone mode and opens a gap in the 
XY fluctuations. By Holstein-Primakoff transformation on the pseudo-spin hamiltonian, 
extending the standard treatment [15] to include d^ff'- 

^0-^^l^|{Wf){deff) 

= 16 4(Jp + Je/2)My|ciC_i| (3) 

where ax = A^x/(2 \M\) is the ordered transverse component of the cr = 1/2 spins (0 < 
Ox < CF — 1/2). The prediction from Eq. (3) apphes to a simple FIMO. In the case of an 
extended FIMO, one expects additional contributions resulting from the mixing of states. 
An "amplification" factor is therefore to be expected. As shown by the dashed line in Fig. 
4, with a multiplying factor ~ 5, Eq. (3) provides a quite good agreement with the data. 
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This factor is surprisingly large and suggests that the dynamics is much more sensitive to 
the extended nature of the FIMO than the statics. 

At H = Hell a FIMO transition is expected to be of second order. The transition 
here is complicated by the frustration: the magnons that soften are incommensurate and 
degenerate in momentum space [5] . The initial ordering vector appears away in wave vector 
space from the minimum gap, which, remarkably, remains non-zero at the transition. The 
complexity seen in our results as to the evolution of the ordering vector towards a locked- 
in commensurate value with coUinear spins, is, we believe, due to interplay of quantum 
fluctuations, enhanced by the frustration of the hexagonal lattice and the anisotropy. Even 
in the commensurate phase a small anisotropy, since it opens a gap, changes the nature of 
the magnon condensation. 

We have observed new and peculiar features at the onset of field-induced magnetic order 
in Cs3Cr2Br9. They should stimulate new developments, concerning simple versus extended 
FIMO, incommensurability and hysteretic behavior. 

We thank E. Ressouche and J. Rodriguez Carvajal for help in the structural refinement. 
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